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Abstract 

An approach to find the field equation solution of the Randall-Sundrum 
model with the IZ2 extra axis is presented. We closely examine the in- 
frared singularity. The vacuum is set by the 5 dimensional Higgs field. Both 
the domain-wall and the anti-domain-wall naturally appear, at the ends of 
the extra compact axis, by taking a new infrared regularization. The sta- 
bility is guaranteed from the outset by the kink boundary condition. A 
continuous (infrared-)regularized solution, which is a truncated Fourier se- 
ries of a discontinuous solution, is utilized. The ultraviolet-infrared relation 
appears in the regularized solution. 
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1 Introduction 



As an approach to explain the mass hierarchy problem, the Randall-Sundrum 
(RS) model[Q, 0] have been taking people's attention both from the 
phenomenologyQ and from the theory]^, |^. The model has, in fact, some 
advantages compared with other approaches such as the Kaluza-Klein com- 
pactification |^, ^ and the (standard) renormalization group approach . 
The most characteristic point is its exponential damping factor (warp fac- 
tor) which could have the possibility of naturally explaining the broadly- 
spreading mass hierarchy ranging from the cosmological constant (10~^^GeV), 
through the weak physics (lO'^GeV), to the Planck mass (lO^^GeV). Fur- 
thermore the recent progress in the AdS/CFT correspondence g |, 0] indi- 
cates the RS-model solution, which is a classical solution in the 5 dim AdS 
space-time, could be regarded as the renormalization trajectory in the 4 dim 
quantum solution. 

We point out, however, an incomplete aspect in most approaches so far. 
They assume the 5-function or ^-function distribution from the outset as 
a form of the classical solution in order to make a ("infinitely-thin") wall 
configuration. Indeed it gives an easy "tool" to analyse the model in some 
limitted situation. It is, however, obscure from the standpoint of the soliton 
(kink) physics and does miss the important role of the "thickness" in the 
regularization standpoint. The configuration, considered in the RS-model, 
generally has a domain wall structure with some finite thickness which is 
determined by the vacua in the asymptotic regions (or the boundary condi- 
tions) and some parameters in the system. In some limitted configuration, 
the thickness approaches zero and the (5-function (or ^-function) appears as 
a well-regularized object. In such a way, we can understand the real mean- 
ing of the limit from the vacuum structure or the system parameters. This 
looks very important especially to understand the problem of the cosmologi- 
cal constant, which is the vacuum energy of the space-time. Needless to say, 
the system configuration should be derived by solving the field equation in 
a proper way. In some reference the thickness was introduced just by 
smearing the assumed (5-function. Such approach loses the real role of the 
thickness. 

Motivated by the above things, a solution of the RS-model, for the one- 
wall case, has been presented jl^ , [T^ . The points are 1) the wall configu- 
ration is obtained as a kink solution of the classical field equation of the 5 
dim AdS gravity; 2) the (5-function limit is specified by some parameters; 3) 
the vacua (asymptotic states), which are necessary to specify the kink, are 
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Fig.l The scalar field for one-wall configuration |12|. y £ R = (—00,00 



introduced by the 5 dim bulk Higgs potential; 4) the stability of the solution 
is guaranteed by the boundary condition. The obtained solution correctly 
gives the basic physical outputs such as the 4D Planck mass {Mpi ~ yjM^jk) 
and the 4D cosmological constant (A4C; ~ —M^k) in terms of the 5D Planck 
mass (M) and the thickness parameter {k). These properties remain valid 
in the present case because we modify (with much care) only the infrared 
boundary condition. The configuration of a wall and an anti-wall is taken 
in the original work[^, but it is now considered unstable. To cure it, one 
approach is to take into account the radion field|14|. Here we point out 
another possibility of creating the wall-anti-wall configuration by taking a 
new infrared treatment which mimics the lattice situation. 

It is well-known that the massless chiral fermion appears as a zero mode 
bound to the domain wall. Now we recall the similar situation takes place 
in the lattice field theory using the 5D lattice[^, In that case, besides 
the fact that the formulation is discrete, some essential differences are there. 
That is, the extra axis y is regularized to be finite —L < y < +L and 
the boundary condtion for the extra axis is taken to be periodic: y 
y + 2L. In this case a wall appears at y = and an anti-wall appears 
at y = L. For every zero mode at y = 0, there is a zero mode of the 
opposite chirality at y = L. Every mode is chirally paired, which shows the 
vector-like nature of the 5D theory. The chiral anomaly in the 4D theory 
is now understood as the flow of the current through the fifth space [p!7|]. In 
the lattice numerical simulation, this configuration is basically taken, with 



some improvements |18, M, pO], and several physical quantities of QCD 



such as the pion mass are numerically calculated 1 21]. In order to realize the 
similar situation in the Randall-Sundrum model, we need the wall-anti-wall 
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configuration. We present a way to make the configuration from the one-wall 
solution of Ref.||. 

A focus here is to clarify some controversial point, that is, whether IZ2 
compactification of the kink configuration(Fig.l) is compatible with the wall- 
anti-wall configuration. On the one hand, property requires the solution 
to be periodic with some finite periodicity in an extra axis. On the other, 
naive expectation implies its behavior, in an asymptotic region (A), does 
not continuously connect with that in (B) of their adjacent period, as far 
as (A) and (B) are different vacua which is required for the soliton (kink) 
configuration. Clearly this is related to the boundary and stability problems. 
The stability is guaranteed by the kink property: the two vacua (A) and 
(B) are related by the discrete (discontinuous) symmetry, $ <-> — $. To 
solve these problems, some close infrared treatment is necessary. Note that 
physical behaviors in both regions (A) and (B) are the same. (Both have the 
same 5D scalar Riemann curvature.) Further note that we do not consider 
a new solution of the kink-anti-kink type. We will first choose a correct 
coordinate where the wall-anti-wall configuration should appear and then 
change the infrared boundary condition for the one-wall solution previously 
obtained. We will show the new boundary condition gives us the other 
(anti-) wall. The present claim is that the (stable) wall-anti-wall solution 
exists by taking the correct coordinate and the new infrared regularization 
proposed here. 

The present model of 5D gravity-scalar system is introduced in Sec. 2. In 
Sec. 3 we change the extra coordinate from the infinite one y G R = (— cxd, 00) 
to a compact one z S (— ^rc, +\rc) in order to obtain the "size" of the extra 
space Tc. In Sec. 4, by imposing the periodic boundary condition, we extend 
the coordinate region of z to R. The new infrared regularization is explained 
in Sec. 5, where the Fourier expansion of (continuous and discontinuous) 
periodic functions is exploited. Truncation of the infinitely expanded terms 
to the finite ones is the key of the present regularization. In Sec. 6, we present 
the numerical solutions of the present model, and using this result, we fix 
all Fourier expansion coefficients numerically. We conclude and discuss in 
Sec. 7. Some appendices are in order to supplement the text. In App.A, 
a compact coordinate, which is different from the one taken in the text, 
is examined. Some advantageous points are noticed. Numerical results of 
Sec. 6 are explained in App.B. They consist of two standard method of the 
numerical calculus: the Runge-Kutta method and the Least Square method. 
In App.C, some simple function, which imitates the solution of the present 
model, is examined in order to clarify the characteristic properties of the 
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Fig.2 The Higgs Potential V{^), (||). Horizontal axis: From (^), 
V{0) = + A > 0, V{vo) = A < 0. 



We take the following 5D gravitational theory with 5D Higgs potential. 



S[Gab, <^] = J <fXy/^{--M''R - -G^^'dA^dB^ - V{^)) , 



where X {A = 0, 1, 2, 3, 4) is the 5D coordinates and we also use the notation 
(X^) = {x'^, y), fJ. = 0, 1, 2, 3. The coordinate X^ = y is the extra axis which 
is taken to be a space coordinate. $ is a 5D scalar field, G = detG^B, R 
is the 5D Riemannian scalar curvature. M(> 0) is the 5D Planck mass and 
is regarded as the fundamental scale of this dimensional reduction scenario. 
V{^) is the Higgs potential and serves for preparing the (classical) vacuum 
in 5D world. The three parameters A, vq and A in V{^) are called here 
vacuum parameters. A(> 0) is a coupling, vq{> 0) is the Higgs field vacuum 
expectation value, and A is the 5D cosmological constant. See Fig.2. It is 
later shown that the sign of A must be negative for the present domain wall 
configuration. The Einstein equation is given by 



Fourier expansion coefficients. 



2 Model Set-Up 
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Following Callan and Harvey|17|, we consider the case that $ depends only 



on the extra coordinate y: $ = ^{y). Because M-dependence can be ab- 
sorbed by a simple scaling (<I> = M^/'^^,vo = M^/'^vq,\ = M~^A, A = 
M^A, = M~^X^,) we may, for simphcity, take 

M = 1 . (3) 

We explicitly write M only when it is necessary. 

3 Infinite Extra Axis and Its Compactification 

We start with the following 5D metric|^, [l^ . 

ds^ = e-^''''y^r]^^dxf'dx'' + dy^ , y G R = (-oo, +oo) (4) 

where r/^,y = diag(— 1, 1, 1, 1). In this choice, the 4D Poincare invariance is 
preserved. The Weyl factor e^^'^^^^ is called "warp factor" and is determined 
by the 5D Einstein equation. The extra axis taken here is an infinite real line 
R = (— cxD, +cxd). The coordinates (X^) = (x^, y) give one wall configuration 
by taking the boundary condition: $(y) — > ±vo , y ±00 and there exists 
a family of exact solution [^]. Now let us move from the y-coordinate to 
another one ^(See Fig.3). 

^ = - = -tanh(My) , -e (--,-) , (5) 

where a free parameter Tc is introduced as the compactification size. (An- 
other compactification is examined in Appendix A.) In the "wall region" 
|y| <^ 1/M (or 1 2; I <C rc), z-coordinate and y-coordinate are almost same 
except a simple factor : z k, y x In the asymptotic (infrared) regions 

|y| > 1/M (or \z\ rc/2), they differ significantly : z/rc ^ ±(i - e~^*^'^') 
as My — > ±00. Without confusion, we may take 

rc = l ■ (6) 

(When Tc-dependence is required, it is easily obtained by the substitution 
z z/rc- ) In terms of the new coordinate z, the line element (Q) reduces 
to 

ds^ = e-^"^-'^r^^^dx^dx'' + 4^2)2 ^^^ ' 

dz = \{l-Az^)dy , ^^{-\\) ■ (7) 
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Fig. 3 Relation between two coordinates, z (compact) and y (noncompact), 
(^). Vertical axis: ^; Horizontal axis: My. 





V 




Phi (z) 


-1/2 










1/2 






-V 






Fig. 4 The Higgs field with the boundary condition @. Vertical axis: ^(z); 

Horizontal axis: z. 



z = ±^ are the points of the coordinate singularity . The 5D Riemann scalar 
curvature is given by 

R = F{-2Fa" + 5Fa'^ -2F'a') , F = 1 - 4z'^ , a' = ^ , (8) 

and it turns out that there are no curvature singularities anywhere, for the 
solution we will consider. Let us consider the case the 5D Higgs field ^{z) 
has the following boundary condition(Fig.4). 

lim <l>(z) ±vo , f > . (9) 

^^±(i-o) 

±vo is the vacuum expectation value in the asymptotic region z ib(i — 0). 
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Fig. 5 The periodically extended Higgs field ^{z). Vertical axis: $(z); 

Horizontal axis: z. 



4 Extra Axis 

In this section we move to the case where the extra axis is through the 
following procedure. If we can properly regularize the coordinate singularity 
at z = ±1 in (see the next section), the coordinate region — | < z < | 
can be extended to R = (— cxd,oo) as follows. 

1. We require the periodic boundary condition : 

$(z) = $(z + l) , a{z)=a{z + l) (or a' (z) = a' {z + I)) ; (W) 

2. Values of <^ at z = ^ + Z are defined as, 

hI + Z)^0; (11) 

3. The universal covering space is taken to be the real number space : 

[-^4^''^ = ^"°°'°°^=^' ^^^^ 

where Z = {0, ±1, ±2, • • •}. See Fig. 5 for the schematic behavior of ^{z). 
Note that we have here newly defined the values of ^{z) at the singular 
points z = ^ + Z (coordinate singularity, not the curvature singularity). 
The points correspond to y = ±00 of the original coordinate, y. ^{y = ±00) 
are not defined in Sec. 3. In ( pT[) we have specified the present treatment of 



8 



y = iboo, and which should be regarded as a part of the present infrared 
regularization. ^ This specification turns out to be important in Sec. 5. 

We furthermore note that the translation invariance y ^ y + c in (^) 
reduces to the periodicity invariance (discrete version of the translation) 
z z + 1. (The situation is the same as the lattice regularization of the 
continuum space.) The lost of the translation freedom c is traded with the 
freedom of the coordinate choice Vc- 

^{z) defined above has the properties: 

PI Piecewise continuous. {Discontinuous at z = ^ + Z.) 
P2 Piecewise smooth, {^'{z) is piecewise continuous.) 
P3 Periodic (5* "^-symmetry ) : z — > z + 1. 
P4 Odd function of z (Z2-symmetry) : ^{z) = — $(— z). 

We call these "<l>-properties" . 

At present, except for the coordinate-singularity points z = ^ + Z, the 
metric (^) is defined for z G R = (—00,00): 

' e-2-(-)r/^., ' 

F{zf 

where F{z) = 1 — Az^ at this stage (soon redefined) and M = {fJ-,z). The 
Einstein equation @ reduces to the following two coupled differential equa- 
tions for ^{z) and cr{z). 

-yFizf = -lFizf^'' + V , 

^ ^ {a'F{z)y = ]^" , (14) 



(Gmn) = I ^ ' (13) 



4F{z)' ' 4 

^ _ T 

dz ' dz 



where a' = = In order to make the above equations periodic in z 



: z ^ z + 1, we must replace F{z) = 1 — 4z^ by its "periodic generalization" 
0: 



Fiz) 



1 - 4z2 for -^<z<^ 



1 - 4z'^ for \z\ > \ 
{z = z' + n , -\<z' < \ , n = ±1, ±2, • • •) 



^ This procedure reminds us of the similar one in the case of making the sphere topology 
(compact) from the space (non-compact) by introducing the point of infinity. 

^ When we take another compact coordinate w, defined in (^^, this "periodic gener- 
alization" is not necessary. 
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F(z) 


1 




-1/2 


1/2 




Fig.6 [Above] The function 1 - 4z2(-l/2 < z < 1/2); [Below] The 
periodicahy generahzed function [1 — 42;^], (p^). Horizontal axis: z. 



[1 - 



(15) 



See . We denote the periodically generalization of a function f{x) as [/(x)]. 
From the definition, F{z) = [1 — Az"^] has the following properties. 

Fl Continuous function of z. 

F2 Piecewise smooth {F'{z) is singular at z = ^ + Z). 
F3 Periodic (5^-symmetry) : F{z) = F{z + 1). 
F4 Even function of z (Z2-symmetry) : F[z) = F{—z). 
F5 Positive semi-definiteness : F{z) > 0. 

We call these properties "F-properties". This process of replacing (1 — 4z^) 
by [1 — 4z^] of ( |T5| ) should be regarded as a part of the present (infrared) 
regularization. 
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For some later use, we present here the Fourier expansion of F{z). 



2 4 ^ (-iV+i 

F{z) = - + -Y.—^^^<'^^^^) ■ (16) 
^ 1=1 * 



This expression can be taken as the definition of F{z) instead of (|l5|). 

We can fix the asymptotic form of o"'(z) using (U) as follows. In this 
process of "periodic generalization", the boundary condition (|9|) is also gen- 
eralized to be 

lim <^>{z) =±vo , vo>0 . (17) 

2-»±(i-0)+Z 

(Note <l>(±i-|-Z) = as introduced in ([Tl|).) We call the asymptotic regions 
{z\z — > ±(^ — 0) + Z} IR-regions and another regions {z\z — > ±0 + Z} 
UV-regions. In the IR-regions, ^ from the above equation, therefore 
a'F{z) — > constant from the second eq. of (0). Furthermore, using the first 
equation, we obtain 




lim Y.{z) = ±uj , uj = \ . (18) 

where S(z) = F{z)a'{z). Note that E(z) = 2^ for -\ < z < \. The 
behavior of the "warp" factor field Ti{z) will be shown to be similar to the 
Higgs field ^{z). {^{z) and ^(z) will be parallelly discussed in Sec. 5 and 
6.) In eq. (|l^) , we notice A should be negative : A < 0. From the first eq. of 
(p!^, we know 

\F{zf^'^ = h:{zf + V>^ . (19) 



From the field equations (|lj) and the boundary conditions (17,18), we con- 
clude Ti{z) and $(z) are odd functions of z. Hence we have 

S(z) = and $(z) = at z = . (20) 

(This condition will be used to solve the field equations (|l^) numerically. 



The boundary conditions (|l7|Jl8|) cannot be taken due to the singularity. 
See Appendix B.) Applying this result to ([l9|), we know \vq^ -|- A > 0. 
Combining the previous result, we obtain fTT 



^1-0^ < A < . (21) 
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It says the sign of A should be negative (anti de Sitter) and the absolute value 
has the upper bound: |A| < ^fo^ • Near the singular points, z — > ib(^ — 0), 
the asymptotic behavior of the line element is given by, from (llF 



l + 2z 



1 - 2z 



(l-4z2): 



:dz^ 



UJ 




(22) 



which shows the singular points are horizons . The power-law of the Weyl 
factor ("warp factor") indicates the scaling behavior of the system when 
z ±(i -0). See Sec.7. 



5 Infrared Regularization at 2 = | + Z 

We come to the most important part of the present regularization. Be- 
fore the presentation, we give here a mathematically well-known fact. The 
periodic step function 9{x) defined by 

{1 2ne<x < (2n + l)e 

x = ne , (23) 

-1 {2n + l)e < X < {2n + 2)e 

where n G Z (see Fig. 7). It has the following properties: 

Tl Piecewise continuous. {Discontinuous at x G e x Z.) 

T2 Piecewise smooth. 

T3 Periodic (5"'^-symmetry) : 6{x + 2e) = 6{x). 

T4 Odd function of a; (Z2-symmetry) : 6{x) = —6{—x). 

T5 Symmetric with respect to the axes x = {±\ + 2Z)e (UV-IR symmetry). 

These are similar to <I>-properties (for ^ = \) except for Property Tl (the 
number of discontinuous points doubles) and Property T5 (UV-IR relation). 
We call these properties "^-properties" . The periodic step function 0(x), 
which is discontinuous, has the following Fourier expansion. 

vr 2/ + 1 e 

(Compare the Fourier expansion of the continuous function F{z), (|l6|) . Main 
changes are, 21 in F{z) is replaced by (2/ -|- l)/e = 2(2/ -|- 1) for e = 1/2, 
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and the coefficient (—1)'"^^//'^ is by 1/(2Z + 1). The discontinuous case is less 
convergent series than the continuous case. ) When we regularize (H) by 
the finite (L) sum, 

9L{x) = -j2-^sm{{2l + l)7T^} , (25) 
vr 2/ + 1 e 

then Ol{x) has the fohowing new properties compared with 9{x): 

TLl Continuous everywhere x G R = (—00,00) (see Fig. 5). Especially 
Ol^x G eZ) = 0. 

TL2 Smooth everywhere. 

Other items 3,4 and 5 are the same as e{x): TL3=T3, TL4=T4, TL5=T5. 
We call these properties "^L-properties" . This simple example characteris- 
tically shows that a discontinuous function can be naturally regularized by 
a continuous function by truncating the infinite Fourier series by a finite L 
sum. L is here regarded as an infrared regularization parameter. The con- 
tinuousness is indispensable for a wall-configuration with finite thickness or 
for a well-defined regularization. The meaning of 1/L is the "thickness" of 
the walls or anti- walls of 9l{x) at x = eZ. ^ The thickness here is purely a 
regularization effect. See Fig. 7. 

With the above fact in mind, we propose here a new regularization in 
order to treat the singularity at z = ^ + Z of the solutions in the previous 
section. First we know ^{z) and S(z) = F{z)a' behave like a periodic 6- 
function at some parameters limit ( the infinitely-thin wall limit). f\ Both 
satisfy the <I>-properties in Sec. 4. Imitating (25) with e = ^, we take, as the 
regularized solution of (0), the following forms for ^{z) and ^(z). 

^L{z) = v^-Y,-^^m{{2l + l)27:z} , 



^ In the "wall region" around the origin j^j <C 1, all L+1 terms equally dominate in 
the RHS of (H) : Bl^x) « ^ ^f^^ zITT^^^^ + l)7rf } = -^{L+l)x. Therefore the thickness 
w can be defined as : |(L + 1)-^ = i, hence we have w = i{L+\) • same thing can 

be said about all "wall-regions" around x £ 2eZ and about all "anti- wall regions" around 
X G e(2Z -I- 1). It is well-known that, in these "wall and anti-wall regions" the truncated 
function 9l(x) most deviate from 9{x) because the neglected terms(high-frequency modes) 
begin to equally contribute with low-frequency ones (Gibbs's phenomenon). 

^ The infinitely thin wall limit is given by Xvq^ — *■ -l-oo (for a given uq and an appro- 
priately chosen A). 
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Fig. 7 [Above] The periodic step function 6{x) ((|2^) or (p4|)) and [Below] its 
regularized function 6l=q{x) ((p5[)). Horizontal axes: x/e. 
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vr 2/ + 1 

where c's and d! s are some numbers to be determined appropriately and L 
is the new regularization parameter which should be taken sufficiently large. 
At present, the way to fix the coefficients, c's and d's, relies on a numerical 



method. (See App.B.l for solving (14) numerically, and see App.B.2 for 



fixing the coefficients by the least square method.) As explained in App.C, 
the behavior of {di} and {q} has three "phases": i) 2/ + 1 <C 1/{47twuv), 
ii) 2/ + 1 PS 1 / {Attwuv) 1 and iii) 2/ + 1 » 1/{4ttwuv), where wuv is the 
thickness around the UV regions (z ~ Z). The critical value / = L* is given 
by the vacuum parameters: 

L* 1 



- VAV = , (27) 
Tc Sttwuv 

where \/ Xvq^ is identified as the (5 dim) Higgs mass mn defined by mn'^ = 
^V"{vo). The (length) scale wuv is an important quantity in the mass 
hierarchy problem. |^ These values L* and wuv are independent of the reg- 
ularization parameter L. This point should be compared with the thickness 
appeared in 6l{x). The condition for the dimensional reduction, from 5 dim 
to 4 dim, is given by 

L* ~ „ > 1 . (28) 

In the present regularized solution (pBj), the UV-IR symmetry (i.e., symme- 
try w.r.t the axes z = ±^ + Z) holds. Therefore another width wjn around 
the IR-regions (z ~ | + Z) is the same as wjjv 

wiR = wuv ■ (29) 

wiR is protected against L(regulator) dependence, at least, taking the reg- 
ularized form of ([2^). (This should be compared with the Ol{x) case, where 
WIR = wuv ~ L~^-) 

Both $i and have the 6'L-properties with e = ^. We consider the case 
that L* is large(L ^ L* S> 1), that is, the solutions $(z) and T,{z) are near 
the ^-function. The infinitely-thin wall limit {6 or 5- function distribution) 
corresponds to, in (pq), the following case. 



L > L* ^ 00 , di^l , ci^l for all Vs . (30) 



' This new mass scale L* jr^ corresponds to the parameter k in the original RS model|l|, 
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Taken into account the condition that the 5D classical Einstein equation 
works, that is, the new mass scale L* /tc should be much less than 5D Planck 
scale M, we should have the following relations between parameters. 

- < — < M . (31) 

rc Tc 

We regard the above three parameters, M(fundamental scale), rc(compactification 
size) and L* (wall-thickness parameter), as the fundamental parameters of 
the theory. 

The numerical results of ( p6|) are given in the next section. (See App.B 
for further detail of the calculation.) The solution (2^) is the regularized 



solution of (14), not a true one . (How to improve (2g) perturbatively, in 
order to approach a true solution, is proposed in Sec. 8. In the practical and 
numerical point of view, the solution ( p6|) is sufficiently close to the true 
solution.) It is, however, sufficient to claim the existence of the solution of 
(|l^ that has the wall-anti-wall (kink-anti-kink) configuration. 



6 Final Numerical Result of jZ^ Compactifica- 
tion 

In Fig. 8, 9 and 10, we plot three sample solutions of ( p6| ) corresponding to 
the following three vacua respectively. The configurations approach to the 
^-function in the order of Vac.l, 2 and 3. 

Vacuum 1 A = 20.0(input) , = l.O(input) , A = -1.88855 

{uo = 1.12207) ,L = 19(input); wuv ~ 8.9 x 10~^(L* ~ 4.5). (32) 
Vacuum 2 A = 40.0(input) , = l.O(input) , A = -3.77762 

(w = 1.58695) ,L = 19(input); wuv ~ 6.3 x 10"3(L* ~ 6.3). (33) 
Vacuum 3 A = lOO.O(input) , vq = l.O(input) , A = -9.4440537 

{to = 2.5091903) ,L = 19(input); wuv ~ 4.0 x 10"^(L* ~ 10.). (34) 

Note that the values of the (5D) cosmological term. A, are very finely 
chosen so that the boundary conditions ( p^ ) and (|18|) are satisfied. All digits 
appearing are "significant figures". As the configuration approaches the 6- 
function limit (mjf^ = Xvq^ — > oo), the necessary number of digits increases. 

^ In ( [31] ) the rc-dependence is explicitly written. The eq.(^), which was introduced 
purely for the notational simplicity, should be taken off and, instead, Mrc 3> 1 should be 
considered. 
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Phi 



0.1 0.2 0.3 0.4 0.5 



dl 

Ir 



. 8 
. 6 
. 4 
0.2 



2.5 5 7.5 10 12.5 15 17.5 



Fig.8 Vacuum 1 (A, vo,A) = (20.0, 1.0, -1.88855): [Top] Higgs Potential 
V{^), (HI). The horizontal axis is <I>. [Middle] Numerical result of ^{z) 
(dotted points, < z < 0.4) for ( [I^ ) and its least-square fit solution ^l{z), 
(p6|), {0 < z < 0.5). The horizontal axis is z. 
[Bottom] The coefficients {d; ; / = 0, 1, • • • , L = 19} of <!>l{z). The 

horizontal axis is /. 



17 




-2 



Phi, PhiL 



Phi 
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0.2 0.3 0.4 0.5 



dl 

Ir 
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0.4 
0.2 



2.5 5 7.5 10 12.5 15 17.5 



Fig.9 Vacuum 2 (A, vq, A) = (40.0, 1.0, -3.77762): same as the figure 

caption of Fig.8. 



18 




Phi 



0.1 0.2 0.3 0.4 0.5 



dl 

Ir 



. 8 
0.6 
0.4 
0.2 



2.5 5 7.5 10 12.5 15 17.5 



Fig.lO Vacuum 3 (A, vo, A) = (100.0, 1.0, -9.4440537): same as the figure 

caption of Fig.8. 
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(See also Vac.lw in App.A.) This shows the cosmological constant is, for a 
given A and vq, dynamically derived in the present framework. Note that A 
is directly related with the 4 dim cosmological constant [|12|. ^ In App.A, we 
give another result for Vac.l with a different compactification coordinate. 

The solutions of T,{z) and T,i{z) are similar to ^{z) and ^l{z). In 
Fig.ll, and ^Liz) are plotted for the case of Vacuum 3. For each 

vacuum, is always closer to the 0-function limit than ^{z). 

The following items can be read from the above output data. 

1. As the Higgs potential has deeper valleies, which corresponds to the 
case that the 5D Higgs mass {rriH = V Xvq'^) becomes larger, $ and S 
approach the ^-function. All coefficients d;'s and q's are expected to 



approach 1 (the limit of (30)). 



2. The wavy region, explained in (pi]), is not so clear in Vacuum 1-3, but 
can be seen in Vacuum Iw of App.A. 

7 Properties of the Solution 
7.1 Wall- Anti- Wall Configuration 

In the previous section, the regularized (numerical) solutions for the wall- 
anti-wall configuration are given. /Z2 compactification is just taking the 
segment [0, ^] for the periodic coordinate z with periodicity 1. The derivative 
of ^l{z) of Vacuum 3 (Fig. 10) is plotted in Fig.l2. Here we see the present 
approach surely gives the wall-anti-wall configuration. The walls appear at 
both ends of the extra axis, not at some middle points in the axis. The 



situation is the same as that in the lattice domain wall|18, |2^. We stress 



the points: 1) the anti-wall is realized by the present IR regularization where 
the discontinuity of ^{z) and T,{z) at the singular points (^ -|- Z) is avoided 
by truncating the infinite Fourier series; 2) the stability of the solution is 
guaranteed by the boundary conditions; 3) UV-IR symmetry is realized in 
the present form of the regularized solution, n We need not the radion field 



which was considered , for the stability, in |||] and was developed in |14]. 



* The cosmological constant problem was reviewed in [g^l- One standard approach 
to the dynamical cosmological constant is the radiative correction. Such scenario was 
explicitly done, using the Coleman- Weinberg mechanism, in p3| . 

^ In ref. , the wall-anti-wall configuration is considered in a modified RS-model and 
UV-IR symmetry is suggested. 
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Si, SiL 




0.5 ■' 

1 

0.1 0.2 0.3 0.4 0.5 

cl 
It 

0.8- • 

0.6- 

0.4- 

0.2- 

2.5 5 7.5 10 12.5 15 17.5 

Fig.ll Vacuum 3: [Above] Numerical result of T,{z) (dotted points, 
< z < 0.4) and its least-square fit solution Sl(-z), (P^), (0 < z < 0.5). 
The horizontal axis is z. 
[Below] The coefficients {q ;l = 0,1, ■ ■ ■ , L = 19} of T,l{z). The horizontal 

axis is I. 
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dPhiL by dz 




Fig. 12 Plot of — where ^liz) is given in Fig.lO(Vacuum 3). The 

horizontal axis is z. 



(In the place of the radion field T^(x)(eq.(13) of ||T|), from eq.(0), the fixed 
function of z, 4/(1 — 42^)^, appears in the present scenario. ) 

7.2 Brane Tensions at the Wall and the Anti-Wall 

We make a remark on a comparative aspect between the wall and the 
anti-wall. Let us consider the "0-function limit" (thin wall limit) mnfc = 
\/Xvo'^rc ~ L* » 1. The behaviors of ^{z) and T,{z) = a'F{z) are considered 
to be the periodic step functions shown in Fig. 13b and d. Correspondingly 
those of ^' and S' can be written as 

^' = voY.{5{z-n)-5{z-n-]-)} , 
nez ^ 

^> = ujJ2{6{z-n)-6{z-n-h} ■ (35) 
nez 

See Fig. 13a and c. The "warp" factor cr(z) behaves as in Fig.lSe. These 
behaviors should be compared with (8)-(10) of Ref.[P. The parts of minus- 
delta-function in eq.(|35|) correspond to the anti- walls. Now we evaluate the 
integrand of the 5D action near z « O(Wall) and z ^ ±^( Anti- Wall). 
Noting (9^<I> = 0, we obtain 

S = J d'xdz^-^^^{-l{-2F^' + 5S2) - If{z)V - ^(c^2 - .0^)^ 
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(a) A Phi'(z) 



-1/2 

-1 \ /' 



(b) Phi(z) 




1 
1 




(e) 


si(z) 1 




1 

V 1 


^1/2 1 




1/2 1 




-ll 







1 ]_! 



Fig. 13 Behaviors at the thin-wall limit. Horizontal axis: z. (a) (b) 
(c) (d) (e) o-(z). Especially, in Fig.l3(e), cr ~ for 

z ~ +0 and cr ~ (tj/4) ln(2/(l - 2z)) + const, for z ~ (1/2) - 0. 
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= J d^x dzL 

where we consider the general curved space g^y{x) for the 4D world, 
i) z « 

F~l , |S|~a; , e"^"^ ~ 1 , |$| ~ , 

L^^{{2u:M^--vo'—)5{z) + \k] . (37) 
4 Tc 3 

The presence of 5{z) term |^ shows the 3-brane is located at z = 0. The 
brane tension is given by 

Tuv = '^^M'^ -^vo^-^ , (38) 



' c 



which is given by the value of u; = ^— 2A/3M-^ and vq. Tjjy is positive for 

1„ 2L* 

4^0 — 



the large value of uj, 2ujAP > jVo"^ — 



ii) z i - 

2 

F ~ 2(1 - 2z) , S ~ w , e"^'^ ~ {j — ^y^^"^ ' $ ~ -yo , 
^ ~ ^V^{{-2u;{l - 2zr - ^(1 - 2zr+')6{z - i) + ^A(l - 2z)-i} 



,(39) 



as far as w > l/rc- 

iii) z « -i + 0(or +i + 0) 

F~2(l + 2z) , , e'^'^ ~ (^-^)+'^/2 ^ $ ^ ^ 

^ ~ ^^/^{(-2a;(l + 2zr - ^(1 + 2z)-+i)<5(z + i) + ^A(l + 2z)-i} 

~0 ,(40) 

as far as a; > l/rc- We see ii) and iii) are continuosly connected at the 
Lagrangian density level. Note that singularities at z = it^ are avoided by 

In (|37| ), we have used relations; S{z) = lima^+o(l/v^27ra) exp{— 2^/2a^}, 5(2)^(2) ~ 
5{z)/\^2Tva, where yphia is regarded as the thickness of the wall, Tc/L* . 
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the warp factors. The above result shows the anti-wah, in the thin wall limit, 
does not remain at the lagrangian density level. We may say the tension for 
the anti-wall is zero. 



TiR = 



(41) 



This view is also confirmed, as far as the value of w is large, by the energy 
momentum tensor ^J—GTmn where Tmn is given by the RHS of eq.(|2|). 



-6cr 



-gg,A\F^'^ + Iv) 



1 



21.5^ ■ 

1 



^g^A'^il - 2z)^-^^+H{z - 1) + A(l - 2z)i-5-i} 
+ 2zf-^^+H{z + i) + A(l + 2z)i-5-i} 



^4 
z ~ . 



F 



4 + 



-4(T 



^{-v^H{z) + 8A) 

^V^g{-^4{l-2zr-H{z 2 

2^\/^{-^(l + 2zY-H{z + \)+ A(l + 2zY~^} 



z ~ 
z ~ 

F ^ F^ > 
z ~ . 



+A(1 -2z)^-3} 



z ~ 



1 



z ~ 



4 + . 



(42) 



for the infrared 



(We have assumed u) > 3/rc in the vanishing of \/—G 
cases. ) 

These results are accepted because we have changed only the infrared 
boundary condition of the one wall solution. The appreciable appearance 
of the anti-wall occurs only at the level of the equation of motion, not at 
the effective action level. On the wall the dynamics of the 4D world are 
operating, whereas it is suppressed (by the warp factor) on the anti-wall. 
From the situation in the 5D lattice simulation, it is quite interesting to see 
how the anti-wall works to provide the other chirality partner for the zero 
mode fermions bound on the wall. 



7.3 Renormalizarion Group Flow 

We make a comment from the viewpoint of the renormalization group in 
the spirit of AdS/CFT. The basic standpoint is to regard the 5 dim clas- 
sical solution as a scaling trajectory [^, ^, in the 4 dim quantum field 
theory. In the present case, it occurs near the horizons. As z ^ ('"c/2) — 0, 
the asymptotic metric (^) says the scalar density "operator" of the 5 dim 
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cosmological term is approximated by 



= v^)^ « A(l - ^)-'^-y-g(x) ^ K{z)^-g{x) , (43) 



where the general curved space gfj,u{x) is considered for the 4D world and 
A(z) is a 4D scalar. We can interpret the meaning of wm of (|2^ ) as the 
(infrared) scaling parameter when the "scale" z approaches the singular 
point rc/2. 

. (44) 

From above results, we obtain 

1 



A(L*) w A( 



\ujrc-l 



''rM-giizn^-mil'^-^r. + l^-^-^r. + l . (45) 

The last quantity corresponds to the (infrared) renormalization group (/?-) 
function for |A|. We know, from the result of Sec. 6, —A is positive and is 
sufficiently large, hence we conclude < 0. The quantity |A| is infrared 
asymptotic free, that is, it decreases as L* — > oo. |^ 



8 Discussion and Conclusion 



The regularized solution (2£) cannot become a true one even when we take 
L = oo, because one of its properties TL5=T5: symmetric with respect to 
z = ^, does not match with the solution except the ^-function limit. In 
order to approach a true solution, as done for the one- wall case [|l2|, we 
must generalize the form of solution (p6D by replacing the constants di and 



Q by z-dependent functions di{z) and c/(z) in the following forms. 

= «^,o + ^ + + • • • 

. , A,2 [z'] , A,4 [Z^] , .... 

Q(^)=A,o + ^(^ + ^^^ + --- , (46) 



In the renormalizable model of 2 dim R -gravity, the same property of the cosmo- 
logical constant is known In the branc world context, a similar result is obtained and 
analysed in a different treatment [EGI ETI pq]. 
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where [z"'] is the "periodic generahzation" of 2" ( in the same way as F{z) 
in Sec. 4). (Compare with eq.(30) of fl^.) In the above we take only even 
powers of z in order to keep the odd function property ( P4 or TL4=T4). 
Note that the above generahzation breaks the UV-IR symmetry (T5=TL5). 
Therefore the present solution, (^) with above generalization (^6[), can be 
regarded as the perturbation around the UV-IR symmetry limit. 

As for the key equations (^) and their generalization (^6|), we can un- 
derstand them by a set of general properties. Let /(x) be a real function 
defined on j; G R = (—00, 00). If f{x) satisfy the following properties: 

Gl Piecewise continuous everywhere. 

G2 Piecewise smooth everywhere. 

G3 Periodic with the periodicity of 1 (^^-symmetry) : f{x) = fix + 1). 

G4 Odd function (Z2-symmetry) : f{—x) = —f{x). (Using the item G3, 
an important property : /(Z) = is deduced.) 

G5 Symmetric with respect to x = (UV-IR symmetry). 

then, the general form can be written as 

00 

f{x)=J2aism{{2l + l)27Tx} , (47) 
/=o 

where {ai} are constants. If we replace the infinite sum by the finite sum 
J2i=o fo'^ regularization, the word "piecewise" in the items Gl and 02 
can always be removed. As a "deformation" of (|47|), at the cost of the item 
05 (UV-IR symmetry), we can generalize the constants {ai} to 

aiix) = aifi + ^[x] + —[x] + --- . (48) 

This generalization produces UV^IR (Planck^TeV) asymmetry. 



Since Hofava-Witten's paper 1 29 1, S^/Z2 compactification (Z2 orbifold) 
becomes popular as a dimensional reduction procedure in the string inspired 
unified models. It gives essentially an wall at one end of the extra axis and 
the anti-wall at the other. The present infrared regularization serves as 
realizing this configuration. 
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Appendix A : Stereographic Compactification 

Instead of the compact coordinate z defined by (^) , we can take another one 
w defined as 

\ w \ , ^ , 

My = tan7r— , < — <- , - cxd < My < cxd . (49) 

Tc 2 rc 2 

See Fig. 14. We take M = rc = 1. The line element (^) is rewritten as 



2 



2 



ds^ = e-2-(-)ry^,dx'^(ix" + ^-—rjdw 

IT 

'^y = 7 — 7 — ^^"^ • (50) 

{cos(7r7i;)}^ 

t(; = ib^ are the points of the coordinate singularity. The boundary condition 
is 

lim ^{w) ±vq , vo > . (51) 

The Einstein equation (|2|) reduces to the same form as (^) but with a 
different F. 

3 , da .n„, ,9 1„ / / d*^ , r, 

3 1 d ^da ^ ^ 1 ^d^ ^2 

4 Fi (w) dw dw 4 dw ' 

2 

Fi{w) = -{cos(7rw)}2 . (52) 
vr 
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Fig. 14 Relation between two coordinates, w (compact) and y 
(noncompact), (^). Vertical axis: My; Horizontal axis: ^. 
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Compared with the case using the coordinate z, the above equations are 
straightforward to the generahzation from S (— ^, ^) to w G R = (—00, 00) 
because Fi^w) is periodic w.r.t. w ^ w + 1. The properties of is the 

same as those of F{z) except a shghtly better situation in the point F2: 
smooth everywhere. Ah procedures in the text are vahd for the coordinate 
w just by replacing F{z) by Fi^w). 

In Fig. 15, we give a sample result for Vacuum Iw (A=20.0(input), vq = 
l.O(input), A=-l. 888810641, a;=l. 122143972). It should be compared with 
Vacuum 1 in the text. The shape of the Higgs potential V{^) is almost 
same as that of Vacuum l(Fig.8), but that of ^{z) is much closer to the 
^-function. ( In accordance with this, more digits are required for the ap- 
propriate value of the cosmological constant A. ) The "wavy" behavior 
(App.C) is recognized in the plot of {di}. The different choice of coordinate 
gives the different behaviors such as the sharpness of $ (or the value of 
wuv)- 



Appendix B : Numerical Results in Sec. 6 

The numerical results of Sec. 6, where the regularized solution ( p^ of the 
field equation (IT^) is given, are obtained as follow. 

Appendix B.l : Numerical Solution by Runge-Kutta 
Method 

First we can directly solve the coupled field equation ( [T^ ) using the 
numerical method. In terms of $(2:) and T,(z) = F[z)a'{z), the field equation 
( |l^ can be written as 

§ = ±;j^{6S= + A(4.= -V)= + 4A)4 , 

where X,vq and A are the vacuum parameters. F{z) is (1 — Az'^) in the 
text (Sec. 4), and (2/7r)(cos vrz)^ in App.A. Due to the periodicity and the 
odd function property, we may focus on an region [0, ^) of z and may take 
only the + sign in the first equation above. The above coupled differ- 
ential equation about {^{z),T,{z)) can be numerically solved, for a given 
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Phi, PhiL 
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0.1 0.2 0.3 0.4 0.5 

dl 

Ir 
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. 6- . 
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Fig.15 Vacuum Iw {\,vo,A) = (20.0,1.0,-1.888810641) : [Above] 
Numerical result of ^{w) (dotted points, < w < 0.4) and its best fit 
approximate solution ^l{w), ([2^). Horizontal axis is w. 
[Below] The coefficients {dr,l = 0,1, L = 19} of <^l{w). 




vacuum parameters (A,t;o,A), by the Runge-Kutta method. As for the 
boundary condition we cannot take (|l^) and (|l8|): hm^_^j__g ^{z) — > vq and 

lim^^i.g S(z) = uj = because of the singularity ai z = \. Instead 

of this condition in the IR (or asymptotic) region, we take 

$(z = 0) = S(z = 0) = . (54) 

This is the condition in the UV (or non-asymptotic) region . Note that the 
above one is the necessary condition when we take the boundary condition 
(P) or the odd function property of '^{z) and S(-z) (see sentences around 
(pO|)). We take the fohowing procedure to have a rehable result. 

1. Taking the above boundary condition, we numerically compute (|53| ) 
for the region [0,0.4]. Here we stop the calculation a little before the 
singularity point z = 0.5. 

2. There appear three patterns in the calculational results: a) They di- 
verge before z reaches 0.4; b) They become imaginary before z reaches 
0.4; c) Calculation lasts to z = 0.4 and all values converge to finite 
ones. Which pattern appears depends on the choice of the vacuum 
parameters {X,vo,k). 

3. Among the results of pattern c), we pick up the best one near to the 
"practical" boundary condition: lim^^o.4 ^{z) vq and lim2^o.4 ^(-2) - 
UJ = ^— 2A/3. Generally, for a given (A, fo)(input), the calculational 
value $ or S tends to be suppressed (or get imaginary ) when |A| 
increases, while it tends to increase ( or diverge) when |A| decreases 
The width of the wall ( or the initial slope ) can be controlled by 
Xvq. After trying the calculation, with 20-50 values of A, for one input 
{X,vo), we can always obtain a satisfactory solution. This says that 
the boundary condition ( |54[ ) and the suitable choice of A realize the 
original boundary condition (^). 

4. For the numerically unaccessible region (0.4,0.5), we put the asymp- 
totic values Vq for $ and 2A/3 for E by hand. Note that this 
region is the asymptotic region, therefore it should be dynamically 
simple although numerically difficult due to the infrared singularity. 

Some sample output data are given in the text {F{z) = 1 — 4z'^) for several 
cases: {\,vo,A) 
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=(20.0,1.0,-1.88855)[Fig.8],(40.0,1.0,-3.77762)[Fig.9], (100.0,1.0,-9.4440537) [Fig.10,11]. 
Another one is given in App.A(F(z) = (2/7r)(cos ttz)^) for (A, uq, A)=(20.0,1.0,- 
1.888810641) [Fig. 15] . 

In ref.[^], using the y-coordinate (Q), some quantities {a'{y),^{y), R,etc) 
are obtained both analytically and numerically. 

Appendix B.2 : Least Square Fitting 

We fit the solution obtained in the previous subsection by the proposed 0- 
th order formula (|2^). L should be taken appropriately large. The critical 
value L* , explained in App.C, can be roughly obtained by 

L 



VW . (55) 



Next we fix the coefficients d's and c's of (pq). A standard way is to minimize 
the following quantity (Least square method). 



I{do,di,---,dL)= f mz)-'^L{z))^dz , 61 = , 
Jo 

J{co,ci,---,cl)= f^^\^{z)-^L{z)fdz , 5J = , (56) 
Jo 



where ^{z) and T,{z) are regarded as the exact solution. Solving this equa- 
tion, we obtain 

-lj^= r'Hz)sin{{2l + l)2nz}dz , (57) 

-K Zi + I Jo 

1 /2 

-— ^ = / J:{z) sm{{2l + l)2Tr z}dz , . (58) 

TT 21 + 1 Jo 

The right hand side of the above equations can be numerically evaluated 
using the numerical results of ^{z) and Ti(z) in the previous subsection. 
The samples of {di} are given in Fig. 8-10 and Fig. 15, and {q} are in Fig. 11. 

Appendix C : Trapezium Model Solution 



As a simplified model of the solution of (14) or (|53|), we can take the 
following simplified model. See Fig. 16. 

^/ N \ when 0< z<2w 

\ vq when Zw < z < ^ 
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Phi ( z ) 






/ ' ' 
/ 1 ' 






w 2w 


1/2 



Fig. 16 The trapezium model solution (|59[). The horizontal axis is z, and 

the vertical one is ^(z). 



We can qualitatively understand the behavior of the coefficients { di } 
appearing in the output data in Sec. 6 and App.A. Using the formula (|5^, 
di is obtained as 

l^ sin{(2; + l)47rzi;} , 
^' = 2^ (2/ + l)4vru; + ' ^ = 0'1'---'^ • (^O) 

Generally three "phases" appear depending upon some regions of I. 



i)0< {21 + l)47rt(; < 1 
if) (2/ + 1)4™ 1 
Mz) {21 + 1)4™ > 1 



di^l- p^w^{2l + 1)^ (parabolic) 
wavy region 
di h (constant) (61) 



The critical value of / = L* is given by 

L* , (62) 

8ttw 

which is independent of the regularization parameter L. In Fig. 17 and 18, we 
plot the above result for (L, i(;)=(19,0.15) and {L,w)= {19, 0.05) respectively. 
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Fig.l7 The plot of {di} for the trapezium model (L, w;)=(19,0.15). 
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Fig.18 The plot of {di} for the trapezium model (|6|). (L, w;)=(19,0.05). 
Axes are the same as in Fig. 14. 
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